
REMARKS 

This amendment responds to the office action mailed on February 19, 2003. In 
the office action the Examiner: 

rejected claims 1 , 3-5, 28, 29, 33, 34, 54, 56, 58, and 60-65 under 35 U.S.C. § 
103(a) as being unpatentable over Tinkham in view of Char ct ai (USP 5,157,466); 

• rejected claims 2, 30, and 31 under 35 U.S.C. § 103(a) as being unpatentable over 
Tinkham in view of Char et al. and further in view of Shnirman et al. (Physical Review B 
57, p. 15400, 1998); 

rejected claims 6, 8-10, 35, 39, 40, 41, 53, 55, 57, and 59 under 35 U.S.C. § 
103(a) as being unpatentable over Tinkham in view of Char ct al. and further in view of 
Baechtold et al (USP 3,953,749); 

rejected claims 7, 1 1-18, 36, 37, 42, 43, 45, 46, and 48-50 under 35 U.S.C. § 
103( a) as being unpatentable over Tinkham in view of Char ct al., Baechtold et al. and 
further in view of Shnirman et al\ and 

* rejected claims 1-18 and 28-65 under the judicially created doctrine of 
obviousness-type double patenting as being unpatentable over claims 1-6, 12, 14, 15, and 
18-25 of U.S. Patent No. 6,459,097 in view of Char et al 

After entry of this amendment, the pending claims are claims 1-18 and 28-65. 
Applicant has amended claims 3 and 65 to correct for antecedent basis. No new matter 
has been added. 

CLAIM REJECTION UNDER 35 U.S.C. § 103(A) 
Legal standard 

When rejecting claims under 35 U.S.C. § 103, the PTO bears the burden of 
establishing a prima facie case of obviousness. In re Bell, 26 USPQ2d 1529 (Fed. Cir. 
1993). To establish a prima facie case, the prior art reference, or references when 
combined, must teach or suggest each and every limitation of the claimed invention. 
MPEP § 706.02(j). The teaching or suggestion to make the claimed invention and the 
reasonable expectation of success must both be found in the prior art, not in the 
Applicant's disclosure. In re Vaeck, 20 USPQ2d 1438 (Fed. Cir. 1991). There must be 
some motivation, suggestion, or teaching of the desirability of making the specific 
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combination that was made by the Applicant. /// re Fine, 837 F.2d 1071, 1075 (Fed. Cir. 
1988). 

In the present instance, one relevant inquiry is whether the cited art, either alone 
or in combination, teaches each and every limitation of the rejected claims. Another 
relevant inquiry is whether the prior art provides one of ordinary skill in the art with a 
suggestion or motivation to modify or combine the teachings of the references relied upon 
by the PTO to arrive at the claimed invention. As discussed in detail below, the cited art 
fails to satisfy either of these requirements. Therefore, the 35 U.S.C. § 103 rejections 
should be withdrawn. 

Claims 1, 3-5, 28, 29, 33, 34, 54, 56, 58 and 60-65 

The PTO has rejected claims 1 , 3-5, 28, 29, 33, 34, 54, 56, 58 and 60-65 under 35 
U.S.C. 103(a) as being unpatentable over Tinkham in view of Char et al. (USP 
5,157,466). Applicant respectfully traverses the rejection. 

One limitation found in each independent claim in the present application is a 
clean Josephson junction having unconventional superconducting material on at least one 
side of the junction, hereinafter referred to as an unconventional superconductor clean 
Josephson junction. In the December 3, 2002 response to the September 3, 2002 Office 
Action, the Applicant argued that the bi-epitaxial Josephson junctions described in Char 
et al. are not unconventional superconductor clean Josephson junctions as recited in the 
claimed invention. In the February 19, 2003 Office Action, the PTO dismissed the 
unconventional superconductor clean Josephson junction limitation in the independent 
claims by stating that it would have been obvious to provide the best quality crystal 
structures since this is standard in semiconductor processing. 

Here, the Applicant will describe the current-phase relationship of Josephson 
junctions made from conventional superconducting materials and unconventional 
superconducting materials. Then, the current-phase relationship of clean Josephson 
junctions will be reviewed. Next, the Applicant will discuss why Char et al. and Tinkham 
do not teach or suggest an unconventional superconductor clean Josephson junction and 
why there is no motivation to modify Char et al. to incorporate clean Josephson junctions 
into Char et al. devices. 
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Joscphson junction current-phase relationships 
In general, the current-phase relation of a Josephson junction is described by an 
odd periodic function commonly represented by the Fourier expansion: 
/((p)^ I x sin(<^)+/ : -sin(2<£>)+..., (1) 

where Ii and h represent the critical current of the first and second harmonics 
respectively. In Josephson junctions formed out of conventional superconducting 
materials, the second harmonic term and higher terms are negligible. See, ITichev et aL, 
1999, Physical Review B 60, p. 3096, (hereinafter "H'ichev 1999") which is enclosed as 
Exhibit A. 

The order parameter of a superconducting material determines the properties and 
characteristics of the superconducting material, and hence the current-phase relationship 
of weak junctions formed in the material. Conventional superconducting materials have 
isotropic order parameters. In contrast, unconventional superconducting materials have 
anisotropic order parameters. See Ifichev 1999. A common unconventional 
superconducting material is the d-wave superconductor YBa 2 Cu307_ x (YBCO), which is 
used in both Char et ai and Ifichev 1999. The term "d-wave" indicates the type of 
symmetry of the anisotropic order parameter. 

Due to the anisotropy of the d-wave order parameter, the current phase 
relationship for a Josephson junction in a d-wave superconductor has a temperature 
dependent second harmonic term. See irichev 1999. The current-phase relationship of 
the Josephson junctions described in Ifichev 1999 is: 

- /; -sin(<p)+/" -sin(2^), (2) 

where l[ and I ! J are the critical currents of the first and second harmonics respectively. 

ITichev 1999 established that the realized non-sinusoidal behavior in the current-phase 
relationship of this unconventional superconducting Josephson junction is explained by 
the presence, and in some cases dominance, of the second harmonic term. 

Clean Josephson Junctions 

The greater the influence of the second harmonic in the current phase relationship 
of a Josephson junction, the greater the deviation from conventional 27T periodic 
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sinusoidal behavior. An unconventional superconductor clean Josephson junction is 
defined by a current-phase relationship in which the second harmonic makes a distinct 
contribution to the characteristics of the junction. In terms of Eqn. 2, this is the regime 
where /" > l[ 12 , which causes the equilibrium state to shift from (p=0, in the sinusoidal 

case, to about ±7t/2, creating a doubly degenerate ground state phase difference across the 
junction. In other words, the phase differences of about +7i/2 and -nil have equal energy 
across the unconventional superconductor clean Josephson junction. 

As has been described in previous office action responses, the doubly degenerate 
ground state is an important feature of unconventional superconductor clean Josephson 
junctions. In the present application, the doubly degenerate ground state is necessary in 
part to create appropriate operational parameters of the claimed devices. For this reason, 
the Applicant positively recites the use of a clean Josephson junction and an 
unconventional superconducting material on at least one side of the clean Josephson 
junction in each of the pending independent claims. 

As seen in Il'ichev 1999, the size of the second harmonic is dependent on a 
number of factors including temperature. Therefore, the second harmonic can easily be 
suppressed (e.g., by raising the temperature of the junction). When the second harmonic 
is suppressed, the junction behaves as a conventional Josephson junction. Parameters 
such as temperature, size, transparency, and roughness, all make important contributions 
to this relationship. In particular, it was shown in Il'ichev et al., 1998, Physical Review 
Letters 81, p. 894, (hereinafter "Il'ichev 1998") which is enclosed as Exhibit B, that for a 
Josephson junction having width of 2\im and operating at a temperature 4.2K, the current 
phase relationship deviated from the standard sinusoidal behavior. It was further pointed 
out in Il'ichev 1998 that this was the first observation of non-sinusoidal behavior in a 
Josephson junction. 

In the present invention, the Applicant exploited the dominant second harmonic 
property of unconventional superconductor clean Josephson junctions to form a useful 
quantum computing device known as a qubit. 
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The combination of Tinkham and Char et ai does not teach each and every 
limitation of the claimed invention. 

Char et ai does not teach or suggest clean Josephson Junctions. Based on the art 
accepted definition of a Josephson junction, the voltage phase relationship of a Char et al. 
device will be sinusoidal in nature if it does not include an unconventional 
superconducting clean Josephson junction. A comparison of Fig. 15 of Char et al. with 
Fig. 4 of Lindstrom et ai, 2003, Physical Review Letters 90, 1 17002 (hereinafter 
"Lindstrom") illustrates the difference between conventional Josephson junctions and 
unconventional superconductor clean Josephson junctions used in devices. Although 
Char et al. describes Josephson junctions formed out of YBCO, an unconventional 
superconducting material, the resulting voltage phase properties of the devices illustrate 
conventional behavior, indicating that the second harmonic is suppressed and the 
Josephson junctions are not in the clean regime. In contrast, Fig. 4 of Lindstrom 
examines the same properties of devices that employ unconventional superconductor 
clean Josephson junctions, and the non-sinusoidal results are quite different from Char et 
al. Lindstrom is attached as Exhibit C. 

As discussed above, Il'ichev 1999 demonstrated that second order harmonics in a 
Josephson junction formed in a thin layer of YBCO is temperature dependent. Thus, the 
voltage phase relationship of a Char et al. device should adopt a sinusoidal waveform at 
high temperatures (68K) and adopt a non-sinusoidal waveform at low temperatures 
(4.2K) if the Char et al. device includes a clean Josephson junction. With this in mind, 
comparison of Fig. 3 of Il'ichev and Fig. 15 ofChar^a/. is instructive. Fig. 3 of ILichev 
1999 shows how the non-sinusoidal behavior of a Josephson junction capable of 
exhibiting second harmonic effects is lost as the temperature of the junction is shifted 
from 4.2K to 40K. In complete contrast, Fig. 1 5 of Char et al. shows that the Char et al. 
device "operates properly" (i.e., exhibits 2ir periodic sinusoidal behavior) at temperatures 
ranging from 4.2K to 77K (see Char et al. column 15, lines 35-40). Thus, Char et ai 
describes Josephson junctions for which the second harmonic is suppressed between 4.2K 
and 77K. This means that Char et al. does not teach or suggest clean Josephson 
junctions. 
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Tinkham does not teach or suggest unconventional superconductor clean 
Josephson junctions. Like Char et ai, Tinkham does not disclose the use of 
unconventional superconductor clean Josephson junctions. 

In the February 19, 2003 office action, the PTO stated that Tinkham shows a 
representative small system (i.e., mesoscopic), made up of a small superconducting island 
connected to charge reservoirs and further, a small superconducting island connected to 
two macroscopic superconducting leads. The PTO noted that Tinkham does not detail the 
materials of the island, the leads or the Josephson junctions. 

Applicant respectfully points out that Tinkham is only considering conventional 
superconducting materials in the devices identified by the PTO. In particular, Tinkham 
does not teach or suggest the use of unconventional superconducting materials in the 
disclosed devices nor does Tinkham teach or suggest the use of clean Josephson 
junctions. 

Tinkham investigates the properties of a small superconducting island connected 
to two charge reservoirs via small Josephson junctions (the S-S-S configuration) in the 
Section entitled "Zero Bias Supercurrents with Superconducting Island and Leads" 
(Tinkham, p. 274). Equation 1.38a of Tinkham (p. 275) calculates the supercurrent in 
such devices as: 



where Is is the supercurrent, e is the charge of an electron, h is Planck's constant divided 
by 27T, and 8 is the sum of the phase differences across the two junctions of the device. In 
the last block of text on page 275, Tinkham comments that this supercurrent equation . . 
agrees with the classical result for two identical Josephson junctions in series...". 

The supercurrent relation that Tinkham calculates for S-S-S devices indicates that 
only conventional superconducting materials are being addressed. For example, the 
current phase relationship that Tinkham calculates depends only on the first harmonic of 
the phase difference across the Josephson junctions. If Tinkham was considering 
unconventional superconducting materials in the S-S-S devices identified by the PTO, the 
current phase relationship in Tinkham would include a second harmonic term as in the 
case of Il'ichev 1999. Since this term is not present in the current phase relationship for 
the S-S-S devices identified in Tinkham by the PTO, it is clear that Tinkham only 
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considers devices consisting of conventional superconducting materials that cannot 
support unconventional superconductor clean Josephson junctions. 

There is no motivation to modify Char et al. or Tinkham to form the claimed devices 

In the February 19, 2003 Office Action, the PTO stated that it would have been 
obvious to provide the best quality crystal structures since this is the standard in 
semiconductor processing. The Applicant respectfully submits that the practice of 
providing the best quality crystal structures would not result in the modification of Char 
et al. or Tinkham to include unconventional superconductor clean Josephson junctions. 
The Applicant has two reasons for this assertion. First, even the best bi-epitaxial 
technology available at the time the present application was filed could not have achieved 
the unconventional superconductor clean Josephson junctions recited in the pending 
claims. Second, even if it were possible to modify Char et al. to make the claimed 
junctions, such junctions would have undesirable electrical characteristics. Because of 
these undesirable electrical characteristics their use in the conventional electronic devices 
described in Char et al. would result in unsatisfactory device performance. This 
reasoning is outlined in the following two subsections. Applicant's reasoning is 
supported by the declaration submitted under 37 C.F.R. § 1.132. However, Applicant's 
reasoning stands on its own without the support of the declaration. 

a. Neither Char et al. nor the best quality crystal structures available for bi- 
epitaxial Josephson junction technology at the time of filing of the application were 
sufficiently advanced to prepare an unconventional superconductor clean Josephson 
junction. 

Bi-epitaxial technology involves forming a grain boundary by creating two 
different substrate regions such that a deposited unconventional superconducting material 
adopts different crystal orientations depending on the respective substrate region over 
which it lies. The substrate regions are prepared on a single-crystal substrate and can 
involve deposition of multiple seed layers in combination with buffer layers and etching 
stages to form the distinct substrate regions. Although this method provides more 
flexibility in the creation of grain boundary Josephson junctions than bi-crystal methods, 
the resulting Josephson junctions are not sufficient to take advantage of some properties 
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specific to the unconventional superconducting materials, such as the second harmonic 
described above. 

The bi-epitaxial technology described in Char et al. illustrates a method for 
creating conventional Josephson junctions, that is, unconventional superconductor 
Josephson junctions having a suppressed second harmonic, as described above in detail. 
The Josephson junctions described in Char et al. are not capable of forming the 
unconventional superconductor clean Josephson junction that is recited in the pending 
claims. 

Further, bi-epitaxial grain boundary Josephson junction technology is complex 
and intricate, and was not sufficiently advanced at the time of filing of the Application to 
form unconventional superconductor clean Josephson junctions. Given the difficulties 
with bi-epitaxial technology at the time the present application was filed, one of skill in 
the art would not have been able to modify Char et al. to produce the devices claimed in 
the instant application. As such, the combination of Char et al. and Tinkham simply 
cannot provide a motivation to modify such references in order to make the claimed 
devices. 

b. Unpredictability of the second harmonic in unconventional superconductor 
clean Josephson junctions. 

Even if bi-epitaxial technology could be used to make an unconventional 
superconductor clean Josephson junction, Char et al. provides no motivation for making 
such junctions. As described above, the current-phase relationship of an unconventional 
superconductor clean Josephson junction has a suppressed first harmonic term and a 
significant second harmonic term. Fig. 4 of Il'ichev et aL, 1999, Physical Review B 60, 
p. 3096 illustrates that the presence of a significant contribution from the second 
harmonic term measurably affects the current-phase relationship of the junction. Further, 
at least in the case of a YBCO thin film with asymmetric 45 degree [00 1 ]-tilt grain 
boundaries, this affect is highly temperature dependent and vanishes at temperatures 
approaching 30 K. Thus, the use of clean Josephson junctions in the devices of Char et 
al. would introduce an unpredictable non-sinusoidal current-phase dependence in such 
devices. Since the devices are typically used in applications such as the precise 
measurement of magnetic fields, this unpredictable current-phase dependence is 
undesirable. 
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The unpredictability in the current-phase relationship of such devices, were they 
to incorporate clean Josephson junctions, would come from at least two sources. First, as 
illustrated in Il'ichev 1999 (e.g., Figs. 3 and 4), the second harmonic contribution 
associated with an unconventional superconductor clean Josephson junction is 
temperature dependent. Second, as detailed in IFichev 1999 and in Lindstrom, state of 
the art methods for manufacturing unconventional superconductor clean Josephson 
junctions still have not developed to the point where the strength of the second harmonic 
can be precisely engineered. Lindstrom is attached as Exhibit C. 

In IFichev 1999, bi-crystal YBCO Josephson junction properties were 
investigated. Films having a thickness of 100 nm were fabricated using standard pulsed 
laser deposition on (001) oriented SrTiC>3 bi-crystalline substrates with asymmetric [001] 
tilt misorientation angles 45° ±1°. They were subsequently patterned by argon ion-beam 
etching into 4x4 mm 2 square washer single-junction interferometer structures as 
illustrated in Fig. 1 of Il'ichev 1999. Five such samples were studied in Il'ichev 1999 in 
which the current-phase relation of each of the samples was compared with the current- 
phase relation of conventional Josephson junctions. Sample No. 1 exhibited the most 
anomalous (non-sinusoidal) behavior, samples Nos. 2 and 3 were less anomalous and 
samples Nos. 4 and 5 did not exhibit second order effects whatsoever. See IFichev 1999, 
page 3097, second column, first full paragraph. Thus, even the state of the art methods 
for manufacturing unconventional superconductor clean Josephson junctions such as 
Il'ichev 1999 presently fail to make clean Josephson junctions with consistent second 
order harmonics. 

Lindstrom fabricated SQUIDs from 250 nm thick YBCO-films deposited on 
SrTiCh-bicrystals. The devices were patterned using e-beam lithography and then 
transferred to a carbon mask employing a multiple step process. Finally, the YBCO was 
etched through the mask using ion-milling. Lindstrom reported that the critical current 
varied from sample to sample (Lindstrom, page. 2, column 2, third full paragraph). 
Further, Lindstrom found that the first order and second order harmonics varied by as 
much as tens times between the two junctions in each of the manufactured SQUIDs. 
(Lindstrom, page. 4, column 1, First paragraph). 

The results of IFichev 1999 and Lindstrom show that each unconventional 
superconductor clean Josephson junction manufactured by such methods would have to 
be characterized to determine the magnitude of the first and second harmonics. Such a 
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step is not presently needed in Char et al. and there is simply no motivation to alter Char 
et al. to introduce such a step since Char et al. does not teach or suggest the use of 
devices that make use of the second order harmonics of unconventional superconductor 
clean Josephson junctions. 

For the reasons presented above, there is no motivation to alter Char et al. in order 
to make devices that have unconventional superconductor clean Josephson junctions. 
Likewise, there is no motivation to alter the teachings of Tinkham to incorporate an 
unconventional superconductor clean Josephson junction. As discussed in detail above, 
Chapter 7 of Tinkham discusses the Josephson effect in terms of first harmonics only. 
Incorporation of an unconventional superconductor clean Josephson junction into the 
structures disclosed in Chapter 7 of Tinkham would be inconsistent with the supercurrent 
relationship (first order effects only) that Tinkham gives for such S-S-S devices. 

For the reasons provided above, Applicant believes that the combination of 
Tinkham and Char et al. does not render obvious claims 1, 3-5, 28, 29, 33, 34, 54, 56, 58, 
or 60-65. Therefore, the Applicant respectfully requests that the rejection be withdrawn. 

Claims 60-63 

Claims 60-63 are patentable over the combination of Tinkham and Char et al. for 
the additional reason that neither reference teaches or suggests the form of quantum 
tunneling recited in these claims. Claim 60 recites a supercurrent that tunnels between a 
first ground state and a second ground state. Tinkham and Char et al., or any 
combination of Tinkham and Char et al., do not teach or suggest an unconventional 
superconductor clean Josephson junction between a mesoscopic island and a bank. As 
described above, the unconventional superconductor clean Josephson junctions introduces 
ground states with about a ± 7r/2 phase difference across the junction. Each of these phase 
differences have the same energy, and hence are degenerate, permitting quantum 
tunneling when the elements of the claims are realized. For this reason, the cited 
references fail to disclose a Josephson junction that is configured so that a supercurrent 
proximate to the Josephson junction alternates between a first ground state having a first 
magnetic moment and a second ground state having a second magnetic moment by means 
of quantum tunneling as recited in claim 60. Claims 61-63 depend from claim 60 and 
therefore are not rendered obvious by the combination of Char et al. and Tinkham for the 
same reasons. 
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Claims 2, 30, 31, and 52 

The PTO has rejected claims 2, 30, 31, and 52 under 35 U.S.C. 103(a) as being 
unpatentable over Tinkham in view of Char et aL (USP 5,157,466) and further in view of 
Shnirman et aL 

Tinkham, either alone or in combination with Char et aL, does not teach or 
suggest claims 1, 28, 60 or 64 for the reasons discussed above. Shnirman merely teaches 
a SET capacitively coupled to a Josephson junction qubit. As such, Shnirman et aL does 
not remedy the deficiencies identified in Tinkham and Char et al. . For this reason, the 
combination of Tinkham, Char et aL, and Shnirman et aL does not render claims 1, 28, 60 
or 64 obvious. Since claims 2, 30, 31, and 52 ultimately depend from 1, 28, 60 or 64, the 
cited references do not anticipate these claims either. Accordingly, Applicant respectfully 
requests that the rejection be withdrawn. 

Claims 6, 8-10, 35, 39, 40, 41, 53, 55, 57, and 59 

The PTO has rejected claims 6, 8-10, 35, 39, 40, 41, 53, 55, 57, and 59 under 35 
U.S.C. 103(a) as being unpatentable over Tinkham in view of Char et aL (USP 5,157,466) 
and further in view of Baechtold et aL (USP 3,953,749). Tinkham, either alone or in 
combination with Char et aL, does not anticipate claims 1, 28, 60 or 64 for the reasons 
discussed above. Baechtold et aL merely teaches a binary circuit consisting of a 
series/parallel arrangement of Josephson junctions. As such, Baechtold et aL does not 
remedy the deficiencies identified in the combination of Tinkham and Char et aL For this 
reason, the combination of Tinkham, Char et aL, and Baechtold et aL does not render 
claims 1, 28, 60 or 64 obvious. Since claims 6, 8-10, 35, 39, 40, 41, 53, 55, 57, and 59 
ultimately depend from 1, 28, 60 or 64, the cited references do not anticipate these claims 
either. Accordingly, Applicant respectfully requests that the rejection be withdrawn. 

Claims 7, //, 12-18, 36, 37, 42, 43, 45, 46, and 48-50 

The PTO has rejected claims 7, 11, 12-18, 36, 37, 42, 43, 45, 46, and 48-50 under 
35 U.S.C. 103(a) as being unpatentable over Tinkham in view of Char et aL, Baechtold et 
aL and further in view of Shnirman et aL As previously discussed, Tinkham, either alone 
or in combination with Char et aL, does not anticipate claims 1 , 28, 60 or 64. Baechtold 
et aL and Shnirman et al. do not remedy the deficiencies identified in the combination of 
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Tinkham and Char ct ai For this reason, the cited references, either alone or in any 
combination, do not render claims 1, 28, 60 or 64 obvious. Since claims 7, 11, 12-18, 36, 
37, 42, 43, 45, 46, and 48-50 ultimately depend from claims 1, 28, 60 or 64, the 
references do not anticipate these claims either. Accordingly, Applicant respectfully 
requests that the rejection be withdrawn. 



The PTO has rejected claims 1-18 and 28-65 under the judicially created doctrine 
of obviousness-type double patenting as being unpatentable over claims 1-6, 12, 14, 15, 
and 18-25 of U.S. Patent No. 6,459,097 in view of Char et aL Although the Applicant 
disagrees with the PTO's reasoning, in order to expedite prosecution, the Applicant has 
authorized the enclosed terminal disclaimer. 

CONCLUSION 

In light of the above remarks, the Applicant respectfully requests that the PTO 
reconsider this application with a view towards allowance. The PTO is invited to call the 
undersigned attorney if a telephone call could help resolve any remaining items. 



DOUBLE PATENTING REJECTION 



Respectfully submitted, 



PENNIE & EDMONDS LLP 



By: 




Gary S. Williams 
Reg. No. 31,066 



3300 Hillview Avenue 
Palo Alto, CA 94304 
Telephone: (650)493-4935 
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APPENDIX A 
AMENDED CLAIM 

3. (Amended) The structure of claim 1, wherein the clean Josephson junction comprises a 
grain boundary between the bank and the island. 

65. (Amended) A quantum computer comprising the quantum register of claim 64 and a 
readout device for detecting whether the supercurrent of each clean Josephson junction 
has the first magnetic moment or the second magnetic moment. 
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EXHIBIT A 
IPichev et aL, 1999, Physical Rev iew B 60, 3096 
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PHYSICAL RI YILW B YOl.L'Mi: 60, NUMBLR 5 1 AUGUST 1 999-1 

Anomalous periodicity of the current-phase relationship of grain-boundary Josephson junctions 

in high-T t , superconductors 

L. ILiehev. V. Zakosarenko, R. P. J. Usselsteijn, II E. Hoenig, V. Schultze, and II. -G. Meyer 

Department of Cryoelectronics. Institute for Physical High Technology, P.O. Box 1 00239. D-07702 Jena. Germany 

M. Grajcar and R. Illubina 

Department of Solid State Physics, Comcnius University, Mlynska Dolina F2, $42 15 Bratislava, Slovakia 

(Received 28 January 1999) 

The current-phase relation (CPR) for asymmetric 45° Josephson junctions between two (/-wave supercon- 
ductors has been predicted to exhibit an anomalous periodicity. We have used the single-junction interferom- 
eter to investigate the CPR for these kinds of junctions in YBa : Cu 3 0- t thin films. A remarkable amplitude of 
the ^-periodical component of the CPR has been experimentally found, providing an additional source of 
evidence for the (/-wave symmetry of the pairing state of the cuprates. [SO 163- 1829(99)05629-5] 



A number of experimental results confirm d x i ,.:-wave 
symmetry of the pairing state of high-temperature 
superconductors. 1 An unconventional pairing state requires 
the existence of zeros of the order parameter in certain di- 
rections in momentum space. Thermodynamic and spectro- 
scopic measurements do indeed suggest their existence, but 
by themselves they do not exclude conventional .v-wave pair- 
ing with nodes. 1 Direct evidence for the </-wave pairing state 
is provided by phase-sensitive experiments, which are based 
on the Josephson effect. 2 Quite generally, the current-phase 
relationship (CPR) of a Josephson junction, I(<p) is an odd 
periodic function of (p with a period 2rr.\ Therefore /(cp) 
can be expanded in a Fourier series 

f(if ) = I] sin (p+ / : sin 2 <pH . (1) 

In the tunnel limit we can restrict ourselves to the first two 
terms in Eq. (1). Since the order parameter is bound to the 
crystal lattice, !{<f) of a weak link depends on the orientation 
of the c/-wave electrodes with respect to their boundary. The 
existing phase-sensitive experiments exploit possible sign 
changes of I\ between different geometries." In this work we 
present a phase-sensitive experimental test of the pairing 
state symmetry of cuprates. Namely, in certain geometries, 
the /, term should vanish by symmetry. In such cases, the 
CPR should exhibit an anomalous periodicity. 

Let us analyze the angular dependence of / , 2 in a junction 
between two macroscopically tetragonal c/-wave supercon- 
ductors. As emphasized in Ref. 4. also heavily twinned 
orthorhombic materials such as YBa 2 Cu :i 0 7 v belong to this 
class, if the twin boundaries have odd symmetry. We con- 
sider an ideally flat interface between two superconducting 
electrodes. Let 0 ] (0 2 ) denote the angle between the normal 
to the grain boundary and the a axis in electrode 1 (2), see 
Fig. 1. If we only keep the lowest-order angular harmonics, 
the symmetry of the problem dictates that 4 

cos 2 0, cos 2 0 2 + /, sin 2 0, sin2 0 2 . (2) 

The coefficients I ( J, are functions of the barrier strength, 
temperature 7\ etc. The I 2 term results from higher-order 
tunneling processes and we neglect its weak angular depen- 



dence. It is seen from Eq. (2) that the criterion for the obser- 
vation of an anomalous period of the CPR, /, =0. is realized 
for an asymmetric 45° junction, i.e., a junction with 0, 
-45° and 0 2 = O. 

The I z term is also present in weak links based on con- 
ventional s-wave superconductors but for all known types of 
weak links \! 2 i 1 x \< \ . For instance, for a tunnel junction 
|/ : /7|l<U. For a superconductor normal-metal 
superconductor (SNS) junction, at T=0/' and the 

Fourier expansion of Eq. (1) leads to Iiil^-— 2/5. There- 
fore a possible experimental observation of |/ 2 1 1 1 in an 
asymmetric 45° junction provides direct evidence of t/-wave 
symmetry of the pairing state in the cuprates. 

We have investigated the CPR of YBa : Cu ;i 0 7 - v thin-film 
bicrystals with asymmetric 45° [001]-tilt grain boundaries as 
sketched in Fig. 1, using a single-junction interferometer 
configuration in which the Josephson junction is inserted into 
a superconducting loop with a small inductance L. In a sta- 
tionary state without fluctuations, the phase difference if 
across the junction is controlled by applying an external 
magnetic flux <L> ( , penetrating the loop: <f = if l . — /3f(if-). Here 
<p e = 2 7r<I> t ./<I> 0 ; O 0 =2.07X 10 15 Tnr is the flux quantum; 
/( (p) = /( <p)// 0 is the CPR normalized to the maximal Jo- 
sephson current / 0 , and f3^2rrLI {) i<l> () is the normalized 
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I Ki. 1. Washer-shaped interferometer with one short Josephson 
junction (not in scale). Dimensions are given in the text 
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critical current. In order to obtain the CPR for the complete 
phase range - rr^ip^ir the condition j3< I has to be ful- 
filled, because for [3> 1 the curve tp(</\ ) becomes multival- 
ued. Following Ref. 3, we express the effective inductance of 
the interferometer using the derivative /' with respect to ip 
as L hf( — L[ 1 + \!pf{ip)\. The inductance can be probed by- 
coupling the interferometer to a tank circuit with inductance 
A quality factor C\ and resonance frequency w {) through 
the mutual inductance A/. s External flux in the interferometer 
is produced by a current I dc + I rf - in the tank coil and can be 
expressed as ip e = 2 tt( I dc + /,.,■) M AI> 0 = <p dc + ip rf . where 
M 2 — k z LL T with k a coupling coefficient. Taking into ac- 
count the quasiparticle current in the presence of a voltage V 
across the junction the phase difference is given by the rela- 
tion ip- ip ik .+ ip rJ — f3f{ ip) - 2 ttt( ip) f7 <I> (1 , where r(ip) 
~L!Rj{ip) with Rj{<p) the resistance of the junction. In the 
small-signal limit <p f . t -<\ and in the adiabatic case wHl, 
keeping only the first-order terms, the effective inductance 
L vf -f of the tank circuit-interferometer system is 



L rff L,\ 



L, 



\=L 1 \ 



k 2 pf{<p) ] 
\+Pf'(<p)} 



Thus the phase angle a between the driving current and the 
tank voltage U at the resonance frequency of the tank circuit 
w 0 is 



tan ac( ip) z 



k 2 QPf'(<p) 



(3) 



Using the relation [ 1 4- /3f'{ <p)]d(p- cltp (Jc which is valid for 
iPrf^l and w7<U. one can find the CPR from Eq. (3) by 
numerical integration. 

The advantage of the CPR measurement of an asymmetric 
45° junction with respect to the by-now standard phase- 
sensitive tests of pairing symmetry based on the angular de- 
pendence of/! is twofold. First, it avoids the complications 
of the analysis of experiments caused by the presence of the 
term f s A Second, flux trapped in the interferometer washer 
(see Fig. 1) does not invalidate the conclusions about the 
ratio \ f 2 \ and hence about the pairing symmetry, which is 
not the case in standard phase-sensitive tests of the </-wave 
symmetry 7 of the pairing state. 

The films of 100-nm thickness were fabricated using stan- 
dard pulsed laser deposition on (001) oriented SrTi0 3 bi- 
crystalline substrates with asymmetric [001] tilt misorienta- 
tion angles of 45°±1°. The films were subsequently pat- 
terned by Ar ion-beam etching into 4;<4-mnr square 
washer single-junction interferometer structures (Fig. 1 ). The 
widths of the junctions were 1 2 ^m, The square washer 
holes had a side length of 50 yum. This geometry of the 
interferometer gives A — 80 p\\. The resistance of a similar 
single junction (without interferometer loop) was measured 
directly and R /> 1 il was found. Therefore the condition 
for the adiabatic limit ojt< 1 is satisfied. For measurements 
of a( ip (h ). several tank circuits with inductances 0.2 0.8 
and resonance frequencies 16 35 MHz have been used. The 
unloaded quality factor of the tank circuits 70< 0< 150 has 
been measured at various temperatures The coupling factor 
A" was determined from the period A/ t/[ of (*(lj c ) using 
MM (/( -<b (} . Its value varied between 0.03 and 0.09. The 
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FIG. 2. Left panel: Phase angle between the driving current and 
the output voltage measured for sample No. 1 at different tempera- 
tures as a function of the dc current I dc . The curves are shifted 
along the v axis and the data for r=30 and 40 K are multiplied by 
factor 4 for clarity. From top to bottom, the data correspond to T 
-4.2, 10, 15, 20, 30, and 40 K. The data measured for 36 3 bicrys- 
tals 36°, #2^0) at 7^40 K in the same washer geometry are 
shown for comparison (open circles). Right panel: The same for 
sample No. 3. From top to bottom, the data correspond to 7*= 4.2. 
10, 15, and 20 K. 

amplitude of /,.,■ was set to produce a flux in the interferom- 
eter smaller than 0.1 <f> 0 to ensure the small-signal limit. 

The measurements have been performed in a gas-flow 
cryostat with a five- layer magnetic shielding in the tempera- 
ture range 4.2^T<90 K. The experimental setup was cali- 
brated by measuring interferometers of the same size with 
24° and 36° grain boundaries. We have studied six samples, 
out of which for four samples the 7r-periodic component of 
l{ip) was experimentally observed. At low temperatures for 
two samples (Nos. 1 and 2) the value of / 2 is larger than / , . 
For sample Nos. 3 and 4 I z is approximately 10-20% of 7 { 
and for sample Nos. 5 and 6 7 2 is negligible. As an example 
we plot the phase angle a as a function of the dc current l iU . 
for sample Nos. 1 and 3 (Fig. 2). The behavior of sample No. 
1 at low temperatures is defined by the tt periodic compo- 
nent of I((p). The curves for sample No. 3 are 2 7r-periodic, 
nevertheless for the curve at 7=4.2 K the local minima 
clearly show the presence of a 7r-periodic component. 

In order to determine the CPR we assume that the period 
of a{I dc ) at 7'= 40 K and A/ Jf = 9.6 ^A, corresponds to 
\ip l/c =-2 7r. We take ip dc = 0 at a maximum or minimum of 
a. This is necessary in order to satisfy /(<p = ()) = 0, as re- 
quired by general principles. 3 The experimentally observed 
shift of the first extreme of cx(lj c ) from !j c = 0 (Fig- 2) can 
be due to flux trapped in the interferometer washer. Most 
probably, this flux resides in the long junction originated by 
the grain boundary crossing the washer of the interferometer. 
This long junction does not play an active role because the 
Josephson penetration depth is much smaller than the junc- 
tion length, and external fields produced by I dc are smaller 
than the first critical field. Nevertheless, the long junction 
sets the phase difference for l ( { c -0 at the small junction. 

In Fig. 3, we show the CPR determined from the data in 
Fig. 2. For all curves we have performed a minimal neces- 
sary shift consistent with 7(<p — 0)~0. Thus we have as- 
sumed that at ifj c ~0 a minimum of a{(p dr ) is realized. For 
an interferometer with a conventional v-wave weak link (and 
also for the 36° junction), at y d , — 0 one gets a maximum of 
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FIG. 3. (a) Josephson current through the junction for sample 
No. 1 as a function of the phase difference determined from the 
data in Fig. 2 The scattering of a(<p) values was reduced by fold- 
ing the data back to the interval (0,7r) and taking the average. Here, 
the symmetry a(<p)~ a{ — <p) was assumed, (b) The same for 
sample No. 3. The data for the asymmetric 36° bicrystal at T=40 K 
(open circles) are also shown. 

a(ip dc ). Note that the minimum of a(cp cA .) at ip lic -0 implies 
a paramagnetic response of the interferometer in the limit of 
small applied fields. 

The amplitude of the -rr-periodic component of the CPR 
decreases drastically with increasing temperature, and at T 
— 40 K its contribution is negligible for all samples. The 
temperature dependence of /, and I 2 could be determinated 
with acceptable accuracy for sample No. 1 only. With de- 
creasing 7", |/ 2 | grows monotonically down to T~4.2 K, 
while the / t component exhibits only a weak temperature 
dependence (Fig. 4). 




20 30 
T(K) 

FIG. 4. Temperature dependence of the Fourier expansion coef- 
ficients /[ 2 determined from the experimental data in Fig. 3(a). 
Solid lines are the Fourier expansion coefficients for the numerical 
data in Fin. 5. 
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FIG. 5. I{ip) calculated according to Eq. (64) of Ref. 11 for a 
junction with 0,=45.5°, 0 2 = O, X */ = 1.5, #c=0.5, and r c = 60 K 
at r=40 K for the 36° bicrystal (open circles) was calculated 
with the same parameters except for # 1 = 36°. 

Our experimental results can be understood as follows. It 
is well known that the microstructural properties of the grain 
boundaries, especiaiy 45° boundaries, are defined by their 
faceted nature. Faceting is an intrinsic property of the grain 
boundaries/ 1 7 and, due to d- wave symmetry of the order pa- 
rameter, the properties of the junctions strongly depend on 
the particular distribution of the facets. Small deviations 
from the ideal geometry of the asymmetric 45° junction lead 
to a finite value of /,. Thus for nearly ideal junctions 
|/ : //j 1 at r— ►O. The region T c can be analyzed quite 
generally within the Ginzburg-L andau theory. Let the elec- 
trodes be described by the (macroscopic) order parameters 
A 1 -2 = I A le'* 1 - 2 . Then the phase-dependent part of the energy 
of the junction is E = a[ A , A * + 1 I.e.] + />[( A , A 2 ) 2 + H.c] 

-! where aji, . . . depend weakly on T n> Thus for T 

close to T c we estimate /, ^|A| 2 *( 7\ - T) and / 2 *|A| 4 
~*(T C —T) 2 t leading to \f 2 /l\\<\. With increasing devia- 
tions from ideal geometry decreases. For large 
enough deviations, negligible values of |/ 2 | are expected 
These expectations are qualitatively consistent with the ex- 
perimental data (see also Fig. 4). 

So far, our discussion was based solely on symmetry ar- 
guments. Let us attempt a more quantitative analysis of our 
data now. Two different microscopic pictures of asymmetric 
45° Josephson junctions between c/-wave superconductors 
have been considered in the literature. The first picture as- 
sumes a microscopically tetragonal material and an ideally 
flat interface. 10 ,: Within this picture, only sample No. 1 can 
be analyzed. Sample No. 2 had /„( T- 1.5 K) = 10 : //A. At 
this temperature only the ^-periodic component of f(<p) was 
observed. At higher temperatures / 0 was not measurable 
/(tp) for sample No. 1 calculated according to the model of 
Ref. 1 1 is shown in Fig. 5. The experimental data can be 
fitted within a relatively broad range of barrier heights. How- 
ever, if we require the /(tp) relation of the 36° junction to be 
fitted by the same (or smaller) barrier height as for the 45° 
junction, we conclude the barrier of the 45° junction to be 
rather low. 14 The dependence of I(<p) on T requires a choice 
of 7\ --(A) K in the non-self-consistent theory of Ref. 1 1 . The 
reduction from the bulk 7\-9Q K is probably due to a com- 
bined effect of surface degradation and order-parameter sup- 
pression at the sample surface. The temperature dependence 
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of the ratio of the tt and 2 tt periodic components in I{<p) is 
seen to be in qualitative agreement with experimental data in 
Fig. 3(a). This is explicitly demonstrated in Fig. 4 where we 
compare the experimentally obtained /, 2 with the results of 
the Fourier analysis of the curves in Fig. 5. The divergence 
of A as 7*— +0 is an artifact of the ideal junction geometry 
assumed in Ref. 11. If a finite roughness of the interface is 
taken into account, this divergence is cut off and the experi- 
mental data in Fig. 4 do indeed resemble theoretical predic- 
tions for a rough interface. 12 However, the non-self- 
consistent theory of Ref. 1 1 is unable to explain the 
experimentally observed steep CPR close to the minima of 
the junction energy [see Fig. 3(a)]. 

In a different approach a heavily meandering interface 
with O^-O^x) is assumed. Now, the critical current density 
y.(.v) is a random function with a typical amplitude 
{\j A X )\)^J c • If tne average critical current along the junc- 
tion (j c )<j c * a remarkable Tr-periodic component is present 
in the CPR. The relation |/ : //,| depends on the distribution 
of j\.{x) and can be much larger than one for <j\ )<€j c . . 15,16 
This model qualitative explains the obtained results for all 
samples, however for a quantitative comparison the actual 
microscopic distribution j L (a) should be known. Note that 
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also within the picture of Refs. 15 and 16 the (/-wave sym- 
metry of the pairing state is crucial, otherwise the condition 
(/ t )<^y c * s difficult to satisfy. 

Our present understanding of /(<p) in the asymmetric 45° 
junction is only qualitative. We cannot say whether the re- 
markable amplitude of the 77-periodic component of I{<p) is 
dominated by the microscopically flat regions. L> or due to the 
spatial inhomogeneity of the junction. This issue requires 
further study. 

In conclusion, we have measured the magnetic-field re- 
sponse of a single-junction interferometer based on asym- 
metric 45° grain-boundary junctions in YBa 2 Cu 3 0 7 v thin 
films. A large ^-periodic component of /(tp) has been ex- 
perimentally found, which is in agreement with theoretical 
predictions for cl x 2. v :-wave superconductors. Hence our re- 
sults provide an additional source of evidence for the t/-wave 
symmetry of the pairing state in the cuprates. 
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For various configurations of Josephson junctions incorporating superconductors with unconventional 
order parameter symmetry, such as most high-T r cuprates, deviations from the standard sinusoidal 
current-phase dependence have been predicted. To this point, these deviations have never been observed 
experimentally. We have measured the current-phase relation of high-TV Josephson junctions, namely, 
YBa 2 Cu307-. v thin film bicrystals, comprising symmetric 45° [001] tilt grain boundaries. The current- 
phase relations of all junctions investigated were found to be extremely nonharmonic, in agreement with 
a ^:_ v 2-wave dominated symmetry of the order parameter. [S003 1-9007(98)06674-5] 



PACS numbers: 74.50. + r 

The current-phase relation (CPR) f(<p) is a character- 
istic property of any weak link connecting two supercon- 
ductors. It describes the dependence of the Cooper-pair 
current I p on the phase difference (p of the order parame- 
ters of both superconducting electrodes. In a general form 
it is expressed as 

/, = -1 1, (i) 

I c being the critical current of the weak link. It was 
shown by Josephson [1] that for ideal tunnel junctions 
between conventional superconductors the CPR is sinu- 
soidal, i.e., f{<p) = sin(^). This sinusoidal dependence 
has been confirmed experimentally numerous times for 
standard tunnel junctions between conventional supercon- 
ductors [2]. 

Recently, it has been revealed that the order parameter 
of most high-TV cuprates is unconventional, dominated 
by a d x 2~ v 2 symmetry component [3-5]. Because of the 
sign change of the order parameter associated with this 
symmetry, strong deviations from the standard sinusoidal 
dependence have been predicted for the current-phase 
relations of various configurations of Josephson junctions 
employing such unconventional superconductors [6-9]. 
In particular, nonharmonic and double-periodic current- 
phase relations are expected for junctions oriented nomi- 
nally perpendicular to the (110) direction of one or of 
both electrodes, as well as for junctions for which the 
(110) direction of one of the electrodes is aligned with 
the (100) direction of the other, such as for 45° [001] tilt 
grain boundaries. 

These predictions are highly unusual. Therefore, an ex- 
perimental clarification of the CPR for high- TV junctions 
for which deviations from a standard harmonic behavior 
are expected is desirable. Such experiments will further 



enhance the understanding of the influence of the order 
parameter symmetry on the properties of grain boundaries 
and high-TV Josephson junctions. In addition, they will 
provide valuable information for the design and use 
of Josephson junction-based circuits, of which many 
characteristics directly depend on the CPR. However, to 
our knowledge, such experiments have not been carried 
out. All available data refer to Josephson junctions 
for which nominally sinusoidal current-phase relations 
are expected. The CPR was measured for weak links 
prepared by ion irradiation [10], for step-edge junctions 
[11,12], and also for 24° bicrystal grain boundaries 
[12]. In nearly all of these cases sinusoidal current- 
phase relations were found. Deviations from a sinusoidal 
dependence have been observed only for one step-edge 
junction, measured at 77 K [12]. These deviations can be 
explained by the influence of thermal noise [13]. 

For these reasons we have investigated the CPR of 
YBa2Cu3C>7- v thin film bicrystals with symmetric 45° 
[00 1 ]-tilt grain boundaries, as sketched in Fig. 1(a). For 
these junctions, strong deviations from a sinusoidal CPR 
are anticipated. 

Following a standard approach [14], the CPR was mea- 
sured using a single-junction interferometer configuration 
in which the Josephson junction is part of a superconduct- 
ing loop with a small inductance L, The phase difference 
<p across the junction is controlled by applying an external 
magnetic flux <I> e penetrating the loop: 

<P = <Pe " Pf(<p) + <Pn + 27TW . (2) 

Here, <p c = 27r<J> e /4>o is the external flux normalized to 
the flux quantum * 0 (= 2.07 > 10" 15 Tm 2 ). The con- 
stant (3 = 2ttLI c /<$>{) is the normalized critical current, 
(p n is a term accounting for the effective noise, and m is 
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/ grain boundary 
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tan a(^ dc ) = 



k 2 QPf{<p{<p * c )) 
1 + Pf'tebp^)) " 



(3) 




jc-facet 



grain boundary 
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FIG. 1. (a) Schematic representation of a symmetric 45° 
[00 1 ]-ti It grain boundary junction in a </<:_ v : superconductor. 
The boundary in the superconducting thin film is meandering, 
leading to the occurrence of tt facets, (b) Schematic of the 
measurement setup. The Josephson junction is denoted by JJ, 
and C indicates the capacitance of the tank circuit. The other 
denotations are explained in the text. 



an integer. As capacitive contributions to the loop current 
are insignificant at the measurement frequencies, the junc- 
tion capacitance has been neglected. The quasiparticle 
current is also negligible, for reasons discussed below. 

The superconducting loop is inductively coupled to a 
tank circuit with inductance Lj [see Fig. 1(b)]. This tank 
circuit is driven with a current I T f at a frequency co and 
a dc current /j c . Thus, (p c can be expressed as a sum 
of a dc and an rf component ip t = ip^ c + <p r f. In this 
arrangement, the effective impedance Z e ff (w) of the loop- 
tank circuit combination is a function of <p c . As shown 
by Rifkin and Deaver [14], the CPR can be obtained 
from this dependence, provided that ip T f <SC 1. To obtain 
the CPR for the complete phase range 0 < cp < 2tt, the 
condition f5 < 1 has to be fulfilled in addition. 

To enhance the accuracy of the measurement, we have 
adapted this common approach and retrieved the CPR 
from a measurement of the <p^ c dependence of the phase 
angle a between the drive current / r f and the tank voltage 
V at the resonant frequency of the tank circuit w 0 . As 
described in Ref. [12], at w {) , the a(<pac) dependence is 
related to the derivative of the CPR f'(<p) = df((p)/d(p 
in the following way: 



Here k is the coupling factor between the tank induc- 
tance and the interferometer, A' 2 = M 2 /(LLt), where M 
is the mutual inductance [Fig. 1(b)]. Using Eq. (3), from 
the measured a(<£Me) dependence /'(^(^dc)) is obtained. 
The CPR is restored by integrating f f ((p(<Pd C )) numeri- 
cally, using the d(p(ip^ c )/d(p^ c dependence obtained from 
differentiating Eq. (2) with respect to </?a c . 

The samples investigated consisted of three bicrys- 
talline YBa 2 Cu 3 07-v films with a T c (R = 0) of 88 K. 
The films, with thickness t = 100 nm, were deposited by 
standard pulsed laser deposition on (00 1 )-oriented SrTiC>3 
bicrystalline substrates [15] with symmetric [00 1 ]-tilt mis- 
onentation angles of 45° ± 2° and were subsequently 
patterned by Ar ion-beam etching into 8X8 mm 2 or 
5 > 5 mm 2 square washer single-junction interferome- 
ter structures. The widths of the junctions were b ~ 
2-3 /xm. The washer holes had a side length of 50 /xm, 
leading to L ~ 80 pH. 

To minimize the influence of external noise, the 
samples were measured in superconducting and double 
magnetic shielding at a temperature of 4.2 K. The 
condition f3 < 1 for the investigated interferometers 
was confirmed experimentally from the character of its 
response versus (p dc [12]. 

YBa2Cu307- r grain boundaries with a symmetric [001] 
tilt angle of 45° typically have a normal-state interface- 
resistivity p n > 1 X 10~ 8 ft cm 2 , which we also mea- 
sured for boundaries fabricated under identical conditions 
as the junctions used in the present experiments. In the 
configuration used, this p n corresponds to normal-state 
resistances R„ > 1 17. Accordingly, the relaxation time 
of the interferometer r = L/R is short (r <5C 1/cun), and 
hence the quasiparticle current is negligible [12]. 

For the measurements of <x((p& c ), two tank circuits with 
quality factor Q = 120 and inductance L r = 0.4 zxH, 
<oq = 30 MHz, and L T = 0.73 /iH, <o 0 = 23 MHz, re- 
spectively, were employed. The phase angle was recorded 
as a function of /j c after amplification of the tank voltage 
by a high-impedance amplifier. The coupling coefficient 
k was determined from the period of the a(I& Q ) depen- 
dence. Values of k = 0.072 and k — 0.054 for the re- 
spective tank circuits are obtained. To ensure the validity 
of the small-signal limit, the measurements were carried 
out with (p r f < 0. 15. 

A typical a(/a c ) dependence is shown in Fig. 2. The 
corresponding CPR, depicted in Fig. 3, is clearly deviat- 
ing from the standard sinusoidal behavior. Samples fabri- 
cated on different substrates and measured with both tank 
coils followed closely the same behavior. It is empha- 
sized that the experimental setup employed and the pro- 
cedure followed are identical to those used to measure the 
current-phase relations of step-edge junctions and thin- 
film bicrystals with a symmetric [001] tilt of 24°. For 
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FIG. 2. Phase angle between the driving current and the 
output voltage measured at 4.2 K as a function of the dc current 
/j L , for an YBa^CirtO/-, single junction interferometer circuit 
containing a symmetric 45° [0()l]-tilt grain boundary. 



all of those samples nominally sinusoidal current-phase 
relations were observed, and all deviations of the apparent 
CPR from a sinusoidal one can be attributed to thermal 
noise [12,13]. 

The measured deviations from a sinusoidal de- 
pendence for the current-phase relations of these 
45° bicrystals are startling. It is important to 
note that the effective Josephson penetration depth 
Aj = [3>o/(47T MoOV)A)] I/2 ~ 5 pm is larger than the 
width of the junction b (narrow-junction limit). Here 
A is the London penetration depth. Although several 
mechanisms are known to cause nonsinusoidal current- 
phase relations for narrow junctions fabricated from 
conventional superconductors, all of these mechanisms 
fail to account for the anomalous dependencies presented. 
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FIG. 3. The normalized current through the junction fif(<f) 
as a function of the phase difference (p restored from the 
measured cv(I^) as shown in Fig. 2. For comparison, the 
function p sin(<p) with f3 = 0.8 is plotted as a solid line. 



First, one potential source of such deviations is thermal 
noise. To evaluate its influence we consider a sinusoidal 
CPR and calculate with Eq. (3) the a(</?d C ) dependence, 
assuming a thermally induced Gaussian spread p((p n ). 
With this, the value of tan a(<Pd C ) is given by [12] 

tana(<p dc ) = k 2 Q/3 

cos (fdpc <p„) 



1 + /3 COS (p{(f c < if n 



p(<Pn)d<p n 



(4) 

Using Eq. (4), minor deviations of the CPR from a stan- 
dard sinusoidal behavior are well described quantitatively 
for 24° boundaries measured at 77 IC [12]. However, no 
realistic set of /3 and <p n exists to account for the large 
deviations of the CPR observed for the 45° boundanes. 

Second, for weak links with a high current density, 
current-induced suppression of the order parameter in 
the electrodes close to the weak link can give rise to 
nonharmonic current-phase relations [16]. It is unrealistic 
that this effect is the cause for the deviations presented 
here, as the intragrain critical current density exceeds 
10 7 A/cm 2 at 4.2 K and is therefore much larger than 
the grain boundary (j c ) < 4000 A/cm 2 at the same 
temperature. 

Third, several additional mechanisms, described in [2], 
lead to deviations from a harmonic CPR. In all these 
cases, the slope of f{(p) at <p = 0 is smaller than at <p = 
7T, which is in contrast to our results. Therefore, these 
mechanisms cannot explain the current-phase relations 
observed either. 

On the other hand, as will be pointed out in the 
following, the measured CPR can be accounted for 
by the unconventional order parameter symmetry of 
YBaoCu^Cb-v and by the microstructural properties of 
the grain boundaries, in particular by their faceted nature 
[17,18]. Interestingly, due to the J v :- v :-wave character 
of the order parameter, the faceting has a more significant 
influence on the electronic properties of boundaries with 
a misorientation close to 45° than on boundaries with 
considerably smaller misorientation angles [18]. This 
also concerns the CPR, as the 45° boundaries contain a 
higher density of facets that themselves show anomalous 
behavior. Two kinds of such anomalies, to be considered 
here, are described in the literature. 

First, due to the sign difference of the adjacent lobes of 
the */ r :- v :-wave order parameter, many facets are biased 
with an additional tt phase shift {tt facets) [17-19]. 
These phase shifts give rise to unconventional junction 
properties, such as a spontaneous generation of magnetic 
flux in the grain boundary junction [19-21]. As described 
in Ref. [21], the local phase difference (f(x) along the 
grain boundary (0 < x < b) can be written as 

<p(x) - ((x) + <Kv), (5) 

where £(.v) is a rapidly alternating function accounting 
for the 77 phase shifts and for the spontaneously generated 
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magnetic flux in the junction, and if/(x) is the remaining 
slowly varying phase difference reflecting the magnetic 
flux in the interferometer loop. For a narrow junction 
(b < Aj), t/j is independent of x. With Eq. (5), the time- 
independent sine-Gordon equation, describing the spatial 
dependence of the local phase difference over the junc- 
tion, becomes 

A3 ^¥- 7 f^«n[* + *<*)]. (6) 

C>*- Or) 

The solution of this equation, and thus the pattern 
of self-generated flux, depends on ip. The redistribution 
of this flux by a change of ift is expected to lead to 
remarkable deviations from a harmonic dependence for 
the CPR measured for the entire junction, also if the local 
CPR is nominally sinusoidal. 

Second, it has been proposed [7,8] that the CPR of 
facets formed by the (110) and by the (100) planes 
of the adjacent grains is periodic with it and thus has 
a double periodicity as compared to the standard case. 
Transmission electron microscopy has revealed that 45° 
[001] tilt grain boundaries in YBaiCu^O?-* tend to be 
composed for a considerable part of such facets [22], For 
the whole junction, this leads to an anomalous CPR: 

/ = / ( .| smift + 7 c o sin 2^ , (7) 

by which the observed CPR can be described. 

In summary, the current-phase relations of grain bound- 
ary junctions with a misorientation of 45° were mea- 
sured with a modified Rifkin-Deaver method. The CPRs 
were deduced from measurements of the phase angle be- 
tween the rf drive current and the rf tank voltage. The 
current-phase relations of the Josephson junctions showed 
pronounced deviations from a harmonic behavior, which 
cannot be accounted for by thermal noise or by other stan- 
dard mechanisms, but are attributed to the ^:__ v :-wave 
symmetry of the order parameter and the faceting of the 
grain boundaries. 
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The predominant r/-\vave pairing symmetry in high temperature superconductors allows for a va- 
riety ofcurrent-pha.se relations in Josephson junctions, which is to a certain degree fabrication eon- 
trolled. In this letter we report on direct experimental observations of the effects of a non-sinusoidal 
current-phase dependence in YBCO dc-SQUIEts, which agree with the theoretical descript ion of the 
system. 
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It is well established [1] that the wave function of a 
Cooper pair in most cuprate high-temperature supercon- 
ductors (HTS) has a f/-\vave symmetry. Its qualitative 
distinction from e.g. the anisotropic s-wave case is that 
the order parameter changes sign in certain directions, 
which can be interpreted as an intrinsic difference in the 
superconducting phase between the lobes equal to n. 

The latter leads to a plet hora, of effects, like format ion 
of Andreev bound states at surfaces and interfaces in cer- 
tain crystallographic orientations [2, 3, 4] The current- 
phase dependence Is(d>) in Josephson junctions formed 
by (/(/-junctions, as well as by sd- junctions comprised of 
a cuprate and a conventional superconductor, depends 
both on the spatial orientation of the d-wave order pa- 
rameter with respect to the interface, and on the qual- 
ity of the latter [5, 6, 7, 8, 9]. Time-reversal symmetry 
can also be spontaneously violated and thus spontaneous 
currents generated [10. 11. 12]. Anothei effect can be 
doubling of the Josephson frequency [6, 13, 14]. 

In this letter we report on expei imental observations 
of strong effects of an unconventional current-phase rela- 
tion on the dynamics of two (-/(/-junctions integrated into 
a superconducting interference device (SQUID) configu- 
ration. 

Since I${6) must be a 27r-periodic odd function, it can 
be expanded in a Fourier series. In the most cases only 
the first two harmonics give a significant contribution to 
the current: 



(i) 



In Josephson systems of conventional superconductors 
the second harmonic will usually be negligible [15] but 
in dd-) unctions the second harmonic may dominate 1 If 
// ; > 1 1/2 the equilibrium state* is no longer 0 = 0 but 
becomes double degenerate at (b ~ ± a.rccos(/ f / /2// / ) — > 
tt/2. The system can then spontaneously break time re- 
versal symmetry by choosing either state. Spontaneous 
currents as well as fluxes can be generated in this state. 



The potential will have the shape of a double well and 
there are reasons to believe that it will be possible to 
observe quantum-coherence in this system. The presence 
of a 2 nd harmonic in the current-phase relation(CPR) of 
a ri(/-j unction was confirmed by ITichev et al [8]. 

A non-sinusoidal CPR of the junctions will change the 
dynamics of a dc-SQUID [16]. Regarding the junctions as 
magnetically small, the super cur rent through the SQUID 
in the presence of an external flux <!>,. = <l>o [& x /27r) can 
be written as 



7 S (0, <;> r ) = I ! cX sin <j> - sin(20) + l l v2 sin(tf> + 0 X ) 



/"sin2(r. 



(2) 



The critical current through the SQUID is given by the 
usual expression I c (<t>x) = max^ / s ( 6, O x ). The time- 
averaged voltage over the SQUID in the resistive regime 
is readily obtained in the resist lvelv shunted junction 
(RSJ) approximation. By introducing 6\<p, o., ] — O2 - <P\ 
and applying the same method as in [17] with the nec- 
essary generalizations, we obtain for the average voltage 
over the SQUID 



U~ J = 



Cm + G 2 



+ g 2 r 

2tt J_ k 



do 



h dd 



d 



13) 



Hero G'1.2 »i'<' t li^ normal conductances of the junc- 
tions, and 



71 1 



d/2) - hio + d/2)) 



0, 



(4) 



gives the difference, 6, in phase drops across each junc- 
tion. In deriving (3,4) we have assumed that the induc- 
tance 1 L is equally divided between the SQUID arms. We 
have also neglected the spontaneous magnetic fluxes in 
the r/</-junct ions, due to their small amplitude [11. 18]. 




0 <i> 

\ I] 



FIG. 1: The loults <»f simulations <»| the /, - o, and I" — <p x 
dependence t,»] lL ,|, -S< <>1 II > wit 1. //, -- 1. //. .. n.l. /.'/ =0.2 
and //■/ = (i 1 (aid. units) The dilfeient t niw.s ioinspond 
to biiis (in i tuts in tin; raiiue / - //, to / - 5//, We assume 
L = 0 and - 

Though ( 1) is tnd\ t'xplit ltlv solvable in the limit A — ► 0. 
it always yield- a[-o. -o,] - ^'[o. o., |. This means that 
the usual u\\ ei si< >n symmetry is letained 

The results of nun i* -ri< . 1 1 < all ul.it n >ns based on CJ) and 
(3) are shown m hg 1 The cusps in the critical » urrent 
correspond to the points at which the global maximum 
in (1) switches from one local maximum to another [16]. 
Note the quasi-'!',, ■ 2-penodicitv of the < urrent inclines 
in the V 6 t pictuic lefleitnig the t mi -cut-phase de- 
pendence (1). and their sin It along the <I>,-axis. which 
depends on the sign of the Lias cunent (as it must to 
maintain the central svmmetrv with respect to the ori- 
gin). The shift does not depend on the magnitude of the 
current since we neglect the self-inductance For large 
biases the 'ho-pei i"di< itv is testojed. Indeed, as the bias 
grows, otic set of minima of the washboaid potential. 
U = [h/2* )[-/' cos,;. 4. ( l n cos 'Jo /oj. disappears 
first unless t he hi st harmonu /' is » ji/d/// /ei«t 

We have tabrn ated and studied a laigc number of dc- 
SQUIDs. The samples were tabi i< ated fiotn 250 nm thick 
YBGO- films de] «o^it ed on Si Ti< > -.-bit i a Is. The grain- 
boundary junt timis (GBJs) an - ot t he as\ mmeti n [001]- 
tilt type with the misi »rienl at n »n angle of 45" (li —45° 
GHJ). Foi more mloi mat ion on GBJs see for example 
reference [ 10]. 

The pattern was defined using E-Bcam lithography 
and then t ransfei ret 1 to a carl m m mask employing <» multi- 
step process. Finally, the YBGO is etched through the 
mask using ion-milling This scheme allows us to fabri- 
cate high-quality bicrvstal junctions as narrow as 0.2 /mi. 




as has been reported elsewhere [20]. In the SQUIDs un- 
der investigation the junctions are nominally 2 ///// wide: 
hence the fabi teat ion-induced damage of the junctions is 
small. 

The measurements were done in an EMC-protected en- 
vironment using a magnetically shielded LHe-cryostat. 
However, the magnetic shielding is imperfect, as is evi- 
dent from the fact, that the expected zero-field i espouse 
of our SQUIDs is not exact ly at zero. The measuring elec- 
tronics is carefully filtered and battery-powered whenever 
possible. In order to measure the dependence of t he crit- 
ical current on the applied field we used a voltage dis- 
criminator combined with a sample-and-hold circuit. All 
measurements reported here were performed at 4.2K. 

The SQUID loops are 15 *: 15 /mi 2 . The numerically 
calculated inductance [21] is approximately 25 pH, yield- 
ing the factor 3 = 2;rZJ r /$o between 0.5-2. 

The SQUIDs were largerly non-hysteretic with a re- 
sistance of about 2 Q. The measured critical current 
varies from sample to sample but is in the range of tens 
of microamperes giving a current density of the ordei of 
J r — 10** A/cm 2 . The estimated Josephson penetration 
length \j — ^o/n/^tt/zoJc \/. is approximately 2 //in in 
all junctions, which means that the junctions are mag- 
netically short. This is supported by the quasi-period 
of the pattern in fig 2 being close to the expected value 
0{)/2X[j^ [17] The differential conductance curves do 
not show any trace of a zero bias anomaly (ZBA), as is 
expected for 0° - 45° GBJs. ZBAs has been observed by 
other groups in GBJs with other orientations [2], 

The critical current is plotted as a function of applied 
magnetic field for two SQUIDs in fig. 3. The result is 
m qualitative agreement with theory if we assume that 
the SQUID junctions have different ratios of the 1 st and 
2 nd harmonics of the critical current. This assumption 
is supported by the fairly small modulation depth (it is 
easy to see from equation (2) that I c would go exactly to 
zero in a SQUID with junctions of identical I C 2/h-\)- 

We can fit the data to equation (2). if w f e compensate 
for the residual background magnetic field and assume 
that we have a small excess current (of the order of a 
few fiA) in the junctions. The fitting parameters again 
confirm that there is a large asymmetry between the arms 
of the SQUIDs. Note, that the model does not consider 
the flux penetration into the junctions. 

ddie result for fields of the order of mT is presented in 
(ig.2 which shows the /(-modulation of the SQUID en- 
veloped by an anomalous Fraunhofer-pattern quite sim- 
ilar to what has been reported by other groups [22, 2'A\ 
for 0° 45° GBJs. Note the inversion symmetry of the 
pattern with respect to the origin. That the global max- 
imum is not m the center can be explained in several 
ways; it has been shown for example that this could be 
<lue to the presence of so-called 7r- loops in the junction 
interface [24]. 

Figure 4 shows the V /i-dependence of one of the 
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FIG. 2: Critical current as a function of magnetic field at 
1.2 K. The dashed box indicates the area plotted in fig. 3a. 




Applied Magnetic Field (|iT) 



(b) 



A "A : A' A <i\ 




\/ i -M!J//.j£L ! W! V 

i r 



Applied Magnetic Field (\iT) 



FIG. 3: Critical cm rent its a function of applied magnetic 
field for two different SQFlDs (hat are nominally identical. 
The solid line represents the fitted expression. The fitting 
parameters are as follows: (a) //, — 9 //.A, l l rl — (1.3 //A, l l r [ — 
3.7 //A and lH = 22.7 /<A (b)I t \ - 7.8 pA. l l rl = 3.0 //.A, 
/// = r>.3 /iA and = ■1.3 //A . In both cases the fit has 
been adjusted with respect to the residual background field 
and the excess current of the junctions. 
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FIG. 4 Voltage modulation as a function of applied magnetic 
field for the SQUID whose - B is shown in fig. 3a. The 
pattern is again inversion symmetric. Note the sign change 
at 100 p\ T which we believe is due to a LC- resonance in the 
SQUID loop. 



SQlTDs. The pattern is again field inversion-symmetric. 
The overall structure is the same as in the model depen- 
dence of fig.l. but there is also an additional shift due 
to self-field effects, which depends on the magnitude of 
the bias current and corresponds (at maximum) to a flux 
~ ().14>(). In a beautiful experiment a similar dependence 
was recently observed by Baselmans et al in a Nb-Ag- 
Nb SNS junction where cm rent-injectois were used to 
change the occupation of current-carrying states in the 
normal region [25] A deviation from the model occurs at 
V r = 100 //V where the minima and maxima switch. This 
is probably due to an LC- resonance in the SQUID. Tak- 
ing L = 25 pH, this would require C -- 0.S pF. which 
agrees with our measurements on single junctions 

Remarkably, the observed offset of the V - B- 
charact eristics with respect to the dinci ion of the bias 
current appears to be a much more robust manifestation 
of the presence of a second harmonic of the .Josephson 
current, than the shape of the I c - B curves itself. We 
observed the shift even in SQIJJDs with the smallest junc- 
tions down to 0.5 //in wide, where the deviations from the 
usual sinusoidal CPR were not obvious from the l r - B 
dependence. 

Generally the nature of the transport through a GH.I 
will depend on its transmissivity D. Il'ichev et al. [8| 
have reported values of D as high as 0.3 m symmet- 
ric: (22.5° - 22.5°) (Id junctions as opposed to the usual 
estimate for a GB.J. D ~- 10 ;> - 10"". Since usually 
/// ■■x D, a high-transnussivity GBJ is required in 
order' to observe effects of the second harmonic. An esti- 
mate of the average transmissivity of our junctions would 
be pail/ A ~ 10" 2 [26] assuming /. the mean free path, 
to be equal to 10 inn and a resistivity in the a-b plane p u \ t 
equal to 10" 1 f?cm. This is still too low to explain the 
strong 2 ,ul harmonic we observe. However, it is known 



from, e.g.. TEM-studies [Hi], that the grain-boundary is 
far from uniform: the properties tan significantly vary 
depending on the local properties of the interface, effects 
such as oxygen diffusion out of the GB etc., which are 
difficult to control. It is therefore reasonable to assume 
that there an 1 many parallel transport channels through 
the CJ B [27, 28]. Channels with high transmissivity dom- 
inate the transport, and might have I) ~ 0.1 even though 
the average transmissivity is much lower. This is also 
consistent with the fact that most of our SQUIDs seem 
to be highly asymmetrical which is to be expected if the 
distribution of channels is random. The ratios of // and 
I J 1 can vary as much as ten times between two junctions 
in the same SQUID, even though the fluctuations of the 
total I c from sample to sample are much smaller. It is 
also (dear from general considerations that a high value 
of l l c l excludes a high value of //, since the 2 ml harmonic 
usually dominates if the odd harmonics of the supercur- 
rent are cancelled by symmetry [29]. 

Recent studies of 0° - 45° GBJs have demonstrated 
that the SQUID-dynamios can be altered by the d-wave 
order parameter in YBCO [30]. It is however important 
to point out that our results do not directly relate to e.g. 
tetracrystal 7T-SQUID experiments; the latter crucially 
depend on having one 7r- junction with negative critical 
current, but still only the first harmonic present in I ( \<P)- 
Our SQUIDs have a conventional geometry, but uncon- 
ventional current- phase relations. 

One explanation for the pronounced effects of the 2 nd 
harmonic could be that relatively large sections of the 
interface are highly transparent and have a low degree of 
disorder. This in turn could be related to our fabrication 
scheme which seems to preserve the integrity of the bar- 
rier. This makes feasible their applicability in the quan- 
tum regime and supports our expectations that quantum 
coherence can be observed in this kind of structures. 

To summarize, we have observed very pronounced 2 nd 
harmonic in the current-phase relation of a 'conventional' 
YBCO dc-SQUID with 0° - 45° grain boundary junc- 
tions It has strongly influenced the SQUID dynam- 
ics. All details of the SQUID behavior were explained 
within a simple model of a dd- junction with relatively 
high transparency. We believe that these effects are im- 
portant for better understanding of HTS Josephson junc- 
tion and SQUIDs 
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